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ABSTRACT 

For a mass-selected sample of 66544 galaxies with photometric redshifts (zphot) 
from the Cosmic Evolution Survey (COSMOS), we examine the evolution of star 
formation activity as a function of stellar mass in galaxies. We estimate the 
cosmic star formation rates (SFR) over the range 0.2 < Zphot < 1-2, using the 
rest-frame 2800 A flux (corrected for extinction). We find the mean SFR to be 
a strong function of the galactic stellar mass at any given redshift, with massive 
systems (log(M/M0) > 10.5) contributing less (by a factor of ~ 5) to the total 
star formation rate density (SFRD). 

Combining data from the COSMOS and Gemini Deep Deep Survey (GDDS), 
we extend the SFRD — z relation as a function of stellar mass to 2; ~ 2. For 
massive galaxies, we find a steep increase in the SFRD — z relation to z ~ 2; for 
the less massive systems, the SFRD which also increases from z = to 1, levels 
off at 2; ~ 1. This implies that the massive systems have had their major star 
formation activity at earlier epochs {z > 2) than the lower mass galaxies. 

We study changes in the SFRDs as a function of both redshift and stellar 
mass for galaxies of different spectral types. We find that the slope of the SFRD- 
z relation for different spectral type of galaxies is a strong function of their stellar 
mass. For low and intermediate mass systems, the main contribution to the cos- 
mic SFRD comes from the star-forming galaxies while, for more massive systems, 
the evolved galaxies are the most dominant population. 

Subject headings: galaxies: evolution — galaxies: starburst: galaxies 
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Introduction 



Among the most outstanding issues in observational astronomy today is understanding 
the physics of the star formation process and the parameters responsible for it. The 
cosmic evolution of the star formation rate (SFR), for example, provides clues towards the 
assembly of mass and the development of the Hubble sequence, while parameters which 
govern the star formation activity are directly responsible for the evolution of galaxies. 
Despite significant progress on this subject in recent years, there are still a number of open 
questions concerning the star formation process in galaxies. 

The main problems in studying star formation in galaxies can be divided into two 
general groups: first, uncertainties due to dust extinc t ion and variation in properties of 



differ ent star formation diagnostics (jCram et al. 



1998 



Hopkins et al. 



2001 



Sullivan et al 



2001 1): second, lack of knowledge about the parameters which control star formation 
activity in galaxies (Ellison et al 2007; Li et al 2007; Kaufmann et al 2007; Gomez et al 
2005). There have been major advances in resolving the problems in the fi rst group by 



comparing star formation rate (SFR) est i mates from differe nt diagnostics fICram et al 



1998 



Hopkins et al 



their difference (ISullivan et a. 



measurements (ICalzetti et al 



2003 



Hod 



<:in£ 



2004 



2001 



2007 



Daddi et al 



20071 ). investigating the cause of 



) and exploring the role of dust extinction in SFR 



Sullivan et al. 



200ll ). Due to severe observational 



constraints, however, there has been little progress in addressing problems in the second 
group. 



Cowie et al 



fll996l ) first studied dependence of the specific SFR (SFR per unit mass) 
with redshift. They found the star forming galaxies to have relatively brighter absolute 
K-band luminosities (Mk) at higher redshifts (2; ~ 1) compared to those at lower-z. 
As the Mk is analogous to the integrated stellar mass of galaxies, they conclude that 
higher mass systems to mainly contribute to the star formation rate density (SFRD) at 
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higher redshifts with the main contribution shifting to lower mass systems at lower- 2;, a 
phenomenon they called " downsizing" . This w as confirmed in later studies of nearby galax ies 



Brinchmann et a. 



2005 



Zheng et al 



2004- 



2007 



Heavens et al 



Scoville et al. 



20041) and more distant systems (jjuneau et al. 



m. 



2007b|). Observed changes in "downsizing" 



function of redshift reveals that, the mo st massive ga. 



proceeding later on to smaller systems (jJuneau et al. 



axies form first, with galaxy formation 



2OO5I). 



The correlation between stellar mass and SFR in galaxies as a function of redshift 



galaxy evolution over cosmic time (jPapovich et al 



2006 



of galaxies and 


lence 


Dickinson et al. 


2004) 



study of these relations for different morphological types of galaxies elucidates the origin of 
the Hubble sequence (Abraham et al 2007; Scarlata et al 2007). 

Here, we study the relation between SFR with stellar mass and redshift, using a 
homogeneous, complete and unbiased sample of galaxies covering a wide range in all these 
parameters. We use the extensive multiwavelength data from the Cosmic Evolution Survey- 
COSMOS (Scoville et al 2007a, Mobasher et al 2007; Capak et al 2007a). The combination 
of the depth, area and wavelength coverage in COSMOS is extremely ben eficial for such 



investigations and complements existing deeper bu t smaller area surveys (jJuneau et al. 



2005 



Glazebrook et al 



2005 



Giavalisco et al. 



20041 ). Results from this study are extensively 



used in Scoville et al (2007b) to explore the influence of local density enhancement on star 
formation activity in galaxies. 

Throughout this paper we assume Hq = 70 kms~^ Mpc~^, Q/^ = 0.70 and Qm = 0.30. 
Magnitudes are given in the AB system. 
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Observations, Photometric Redshifts and Stellar Masses 



We use the photometric data obtained from our imaging surveys of the COSMOS field. 
This includes data in U (CFHT), BVriz (Subaru) and (Kitt Peak/CTIO) bands, covering 
the entire 2 sq. deg. area of the CO SMOS field. Details of the photometric observations 
and data reduction are presented in 



Capak et al. 



(12007af ) and Taniguchi et al (2007). For 



the present study, we use a subset of the i-band selected catalog, extended over the entire 
COSMOS field and complete to i < 25 mag. To this magnitude limit, we derive accurate 
photometric redshifts and spectral types (Mobasher et al 2007) and can reliably estimate 
the SFRs. This covers a significantly larger area with a depth similar to previous studies 
(see the compilation by Hopkins and Beacom 2006). 

Photometric redshifts for galaxies in COSMOS are calculated using the template fitting 
method incorporating priors, as described in Mobasher et al. (2007). For each object, we 
derive the photometric redshift, the redshift probability distribution an d the best-fitting 
spect ral type. The spectral types cove r E, Sbc, Scd and Im templates (IColeman et al. 



1980 



exte nded to the UV and NIR by 



Bolzonella et al. 



templates (iKinney et al. 



2000l ) in addition to two starburst 



19961 ). We also consider extinction as an independent variable in 



the SED fits. 

It is possible that by considering extinction as a free parameter in the fit, we increase 
the likelihood of degeneracy in the final result (i.e. redshift, dust extinction and presence 
of the old population all make the SEDs redder). This is examined in Mobasher et al 
(2007) by studying the accuracy of photometric redshifts by comparing the photometric 
(zphot) and spectroscopic (zspec) redshifts for a sample of galaxies in the COSMOS with 
available such data, expressed by the ratio: = (l^^phot — -2spec|/(l + -z^spec))- The smallest 
rms scatter between the estimated photometric and available spectroscopic redshifts is 
obtained when treating the extinction as a free parameter (see Table 4 and Figure 5 in 
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Mobasher et al 2007). We find = 0.03 with a small fraction (~2.5%) of outliers, defined 
as > 0.3. Since Mobasher et al (2007) use a sub-sample of the COSMOS galaxies with 
available spectroscopy to calibrate photometric redshifts, it is possible that the photometric 
redshift sample is biased against fainter galaxies, for which spectroscopic data do not exist 
in COSMOS. To minimize this, we constrain the sample to galaxies with z < 1.10 as these 
have the most reliable photometric redshifts and least number of outliers (Mobasher et al 
2007). 

We estimate the stellar mass corresponding to each galaxy, using the photometric 
redshift, the rest-frame V-band absolute magnitudes (My) and M/Ly ratio dependent on 
the spectral type estimated from the best fit SED for each galaxy, as described in Mobasher 
et al (2007). 



3. Star Formation Rates from UV Luminosities 

The SFRs associated with individual galaxies are estimated using rest-frame 2800 A 
fluxes. We estimate this for each galaxy, using photometric data covering UBViizK bands. 
This is done by interpolating between the two bands that straddle rest-frame 2800 A 
wavelength, making use of the shape of the best-fit SED when doing the interpolation. 
The exception is at low redshifts (0.20 < z < 0.34), where we use the observed U-band 
magnitude, extrapolated to 2800 A using the best-fit SED. Having the observed apparent 
flux, K-correction and photometric redshifts obtained from the best-fit SEDs, we then 
estimate the a bsolute flux in rest-f rame A =2800 A. A detailed description of this technique 



is presented in 



Dahlen et al. 



(120071). 



The UV luminosity at 2800 A is mainly produced by short-lived O and B stars and 
is therefore closely related to the ongoing star formation activity. Here, we follow the 
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approach in iDahlen et al.l (l2007l ) to relate the UV flux to the SFR. This involves using 
predictions from stellar synthesis models and assuming parametric forms for the past 
star formation history to mimic the evolution of the SFR with redshift, as summarized 
below. We derive the conversion factor between t he UV flux and the ongoi ng SFR, using 



the stellar population synthesis code GALEXEV ( iBruzual &: Chariot 



20031). We assume 



solar metallicity and a modified Salpeter I]V 



turnover at low mass (IBaldry &: Glazebrook 



spa nning 0.1 < M/Mq < 100, but having a 



20031 ). To model the dependence of SFR on 



redshift, we use a parametric fi t to the data fro r a the Great Observatories Origins Deep 



Survey (GOODS), presented in 



Giavalisco et al. 



( 120041 ) and parameterized in 



Strolger et al. 



( I2OO4J ). Convolving the past SFR history with the results from the stellar synthesis models, 
we derive a redshift dependent conversion factor, kx{z), between the ongoing SFR and UV 
flux as 



^2800 



^2800^) 



SFR 



—ergs s 



Hz-^ 



:i^ 



Mq yr- 

The dependence of /c280o(^) on redshift is shown in Figure [H Note that we use the 
same volume averaged past SFR history to calcinate the conversion factor for all galaxies. 
Therefore, galaxies that have SFR histories significantly different from the average, may 
be assigned a SFR that differs from the actual value. Furthermore, galaxies of different 
stellar masses are characterized by different star formation histories (Panter et al 2007), 
introducing a mass-dependent bias into the SFR estimates for individual galaxies. We do 
not correct the SFR calibration here for this effect, as this is highly model-dependent and 
introduces more free parameters into the analysis. Moreover, any such correction would 
complicate the interpretation of a mass-SFR-redshift relation discussed in the following 
sections. We show the redshift dependence of the conversion factor in Figure 1, illustrating 
only a weak dependence with a change of only 2% over the redshift interval considered 
here. This suggests that the dependence on the assumed SFR history is weak when 
deriving the SFR. While this is an acknowledged source of uncertainty affecting calibration 
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of the SFR for any individual galaxy, it has no effect on the estimates of the total SFR 
density (SFRD). This is because the SFRD is calculated by integrating over UV luminosity 
functions estimated from the data, after scaling the UV luminosity appropriately to SFR 
using the redshift dependent conversion shown in Figure 1. We will return to the discussion 
of the accuracy of the SFR calibration again in section 5.1, when we compare our estimated 
SFRDs with those derived independently and compiled from literature. 

Accurate estimate of dust extinction is essential in correcting the SFRs measured from 
rest-frame UV (2800 A) flux. Using the extinction values esti mated for invidual ga laxies 



2000|), we 



through the SED flts, and assuming a Calzetti extinction law (iCalzetti et al. 
derive the reddening, Eb-v, for each galaxy. To examine the accuracy of the Eb-v values in 
this study, we explore their distribution for different spectral types of galaxies. As expected, 
we flnd the median extinction (Eb^v) values to increase in galaxies from old and evolved to 
young and star-forming spectral types. Furthermore, we estimate the dust corrected SFRDs 
as presented in section 5.1, with different extinction corrections: the median Eb-v values 
for each spectral type, individual extinction corrections for each galaxy and a constant 
Eb-v value for all the galaxies. We flnd that the flnal results do not depend on the speciflc 
prescription for extinction correction. 



4. Sample Selection and Completeness 

The sample considered in the present study is selected in z-band. This is the fllter with 
the deepest COSMOS photometry and provides a sample most like a mass-selected survey 
(at least to 2; ~ 0.5). Samples of star forming galaxies, selected based on emission lines 
([Oil], Ha), broad-band (i.e. UV, U- band), far- infrared or radio flux, can be very biased, 
since they select only a particular class of galaxies. This is the main reason for selecting a 
sample resembling a stellar mass-selected survey over much of its redshift range. The flnal 
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sample has: 

• Mb < 25: the majority of the objects fainter than this hmit escape detection at some 
wavelengths (especially in the near- infrared), resulting in poor wavelength coverage of 
the SEDs and thus, uncertain photometric redshifts, SFRs and stellar masses. 

• 0.20 < z < 1.10. The lower-limit for the photometric redshift is adopted to minimize 
extrapolation to rest-frame 2800 A for low redshift galaxies, while the upper-limit is 
selected to avoid objects with uncertain photometric redshifts (Mobasher et al 2007). 

• 9.5 < log(M/M0) < 11.5: this allows a sample complete in terms of stellar mass. 
(Figure [2] shows the distribution of stellar mass as a function of redshift). We find that 
galaxies with redder SEDs (i.e. those with spectral types resembling the early-type 
systems), are highly incomplete to log(M/M0) ~ 9.5 while this is the appropriate 
completness limit for galaxies with bluer SEDs (corresponding to later-type systems). 

• Mq/jt < SF R2S00 < 100 Mq /jt: this will exclude sources with excessive (and 
probably uncertain) SFRs, which could dominate the SFR densities, and AGNs 
wrongly classified as star-forming galaxies due to their extreme UV flux. Figure [3] 
shows the SFR distribution as a function of redshift. Very few galaxies are found to 
have SFR2800 > 100 MQyr~^ and therefore, the results are not sensitive to the upper 
limit for the SFRs. 

• My < —19: to allow the selection of galaxies from the same part of the luminosity 
function at different redshifts and minimize luminosity (stellar mass) dependent 
biases. This criterion excludes only an additional 61 galaxies on top of the above 
criteria, predominantly located in the middle two redshift bins. 

After applying the above selection criteria, we have a total of 66544 galaxies in our 
sample. Following Scoville et al (2007b), we divide the sample into the following redshift 
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bins: 0.20 < z < 0.43 (5594 galaxies); 0.43 < z < 0.65 (9608); 0.65 < z < 0.88 (22374) 
and 0.88 < z < 1.10 (28968). In terms of stellar mass, our combined sample is complete to 
M ~ 3 X 10^ Mq over the redshift range 0.20 < ^ < 1.1. 

The fraction of the AGNs in the final sample is < 1% (a total of 1865 sources. M. 
Salvato -private communication), as identified by their X-ray flux. We remove all the 
identified AGNs from the sample. There is a possibility that the final sample is still 
contaminated by a small number of AGNs. However, while massive galaxies at z ~ 2 could 
host obscure AGNs (Daddi et al 2007), the fraction of galaxies at 2; < 1 with strong nuclear 
activity is small. It is also likely that we miss a population of highly dusty star-forming 
galaxies with masses above our stellar mass completeness limit. Our estimated global SFRs 
are therefore likely to be lower limits, although our corrections for dust obscuration to the 
rest-frame UV flux and the technique used to estimate the total SFRDs are expected to 
account for the majority of such missing systems. 



Results 



5.1. Measurement of the Star Formation Densities 



The aim of this study is to explore the effect of stellar mass on the star formation 
activity in galaxies and to follow its behavior as a function of look-back time. Previous 

and SPSS- Gora ez et al 2005) or. 



20051 ). therefore, suffering 



such studies have either been too shallow (i.e. 2d 
covered only small volumes with very few sources (jJuneau et al.l 
from cosmic variance. The present study, based on the COSMOS field, uses data in a wide 
(2 sq. deg) area to medium depth, providing a highly homogeneous sample in the range 
0.20 < z < 1.10, with well-known selection bias. 



We estimate the star formation rate densities (SFRD) in four redshift intervals: 



- 13 - 



0.20 < z < 0.43; 0.43 < z < 0.65; 0.65 < z < 0.88 and 0.88 < z < 1.10. In each redshift 
interval the SFRD is estimated in two ways: (1) by summing over the extinction corrected 
SFRs for galaxies above the stellar mass completeness limit in each redshift interval and 
normalizing them to the volume corresponding to that redshift slice; (2) by integrating the 
rest-frame UV luminosity function of galaxies (corrected for extinction) in the above redshift 
intervals, measuring the UV luminosity densities and converting them to SFRDs. The 
SFRDs from the two methods are listed in Table 1, with the evolution of the SFRD with 
redshift, derived from the COSMOS galaxies, presented in Figure HI T his is compared with 



the S FRDs estimated independently from other studies, as compiled in ( IHopkins fc Beacom 



20061). All the data i n Figu re 4 are scaled to a modified Salpeter IMF described by 



Baldry fc Glazebrookl (120031 ). Errorbars are measured assuming Poisson statistics. The 
SFRDs per redshift bin estimated from the sum of the contributions from individual galaxies 
is lower than those derived from the rest-frame UV LF (Table 1 and Figure 4). This is due 
to the absence of sources fainter than the flux limit of the survey from this method. The 
agreement between the SFRDs estimated here using the extinction corrected rest-frame UV 
LF, and those in the literature over the redshift range 0.20 < z < 1.10 is very good, given 
that the SFRDs are measured from different diagnostics, different surveys and are based on 
different prescriptions for extinction correction. The agreement here lends support to our 
estimated values for redshift, SFR and extinction. Because of their selection wavelengths, 
however, all these s urveys are biased against very dusty, extremely star-forming galaxies 



( lAfonso et al. 



20031 ). implying the estimated SFRDs in Figure H] are likely lower limits. We 



further discuss this and its implications in section 6. 
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5.2. Relation between the SFRDs and Stellar Mass 

The main result of this paper is presented in Figure [5], where we find changes 
in the SFRDs as a function of redshift and stellar mass. We divide the sample into 
three mass intervals: 9.5 < log(M/M0) < 10.0, 10.0 < log(M/M0) < 10.5 and 
10.5 < log(M/M0) < 11.5 in order to have a sufficient number of galaxies in each stellar 
mass interval. The points corresponding to the lower mass bins are likely to be lower limits 
due to incompletness. Errorbars correspond to Poisson statistics, and are often smaller than 
the symbol size. We find that for galaxies in all mass intervals, there is a clear increase in the 
SFRDs with redshift, changing by a factor of 2 to 4 over the redshift range considered here. 
Moreover, at any given redshift, more massive galaxies (log(M/M0) > 10.5) contribute less 
(by a factor of 5) to the total SFRD. This is because they have already gone through their 
intensive star formation phase and have built up their mass. Results from Figure [5] also 
imply that, while massive galaxies continue to undergo modest star formation activity to 
2; ~ 1, they have acquired most of their mass before this redshift (i.e. z > 1). The results 
here are consistent with the "downsizing" pic ture of galaxy formation, where the SF R 



changes from high-mass to low-mass systems ( jCowie et al. 



(119961); 



Juneau et al. 



(120051)). 



The mass-dependent evolution of the SFRDs show progressively smaller contributions 
to the total SFRD (including all masses) as the SFRD declines to \ow-z. The fraction of 
the total SFRD, however, is more or less constant over the redshift range probed for each 
mass bin (i.e. at any given mass bin, the ratio of the SFRD in that bin to the total SFRD 
does not change over the redshift range considered here). This result differs significantly 
from that found by Seymour et al (2008) for the evolution in the high-luminosity end of the 
luminosity function, where very strong mass-dependent evolution is seen for this fractional 
contribution to the total SFRD. This result is likely to be due to the different region of the 
luminosity function sampled by different selection effects arising at different wavelengths. 
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Furthermore, the radio selected sample used in Seymour et al (2008) is mainly dominated 
by galaxies from the high-end of the luminosity function while the optically selected sample 
here is more typically dominated by L* galaxies. 

Using a sample of galaxies detected at 24 /xm wavelength, Zheng et al (2007) studied 
the evolution of the SFRD with redshift as a function of stellar mass. Their sample covers 
the same redshift range as the present study and uses a "dust-free" measure of the SFR by 
combining the UV and total infrared flux (8-1000 jim) for individual galaxies. Therefore, 
by comparing our results with those in Zheng et al (2007), one could quantify possible 
dust-induced biases in Figure 5, despite the very different selection criteria used for these 
samples. The results from the two studies are compared in Figure 6, where they are 
divided into the same mass intervals. Over the redshift range covered by these studies, 
the rate of change of the SFRDs with redshift (i.e. the slope of the SFRD- 2; relations) 
is the same, independent of the stellar mass. For the intermediate mass galaxies with 
2Q10.19 ^ M/Mq < lO^-*^-^^, there is excellent agreement between the the two studies. 
However, for the lowest mass bin at 10^-^ < M/M© < 10^°•^^ the SFRD estimated here 
is significantly higher (by a factor of 5) at all redshifts. This discrepancy is not caused 
by dust extinction in our estimate of the SFRD compared to Zheng et al (2007), as this 
would have led to a relatively smaller SFRD (with respect to Zheng et al) for our sample. 
This is likely a consequence of both incompleteness and selection effects in the lowest 
mass bin of Zheng et al (2007). Their method of infrared (24 //m) stacking, based on the 
optically selected sample of COMBO- 17 sources in this bin, is biased against optically 
faint low-mass galaxies. These are exactly those systems that are likely (by virtue of dust 
obscuration) to be, on average, brighter in the infrared wavelengths than that measured by 
the stacking result from the optically detected sources. This causes a bias, leading to an 
underestimation in the inferred SFRD for their lowest mass bin. Moreover, the inclusion of 
low-mass red galaxies in the stacking analysis (as done by Zheng et al 2007) would serve 



- 16 - 



to reduce the infered average SFR, adding another bias towards underestimating the true 
SFRD. A detailed investigation of the extent of the bias in the low mass bin of Zheng et al 
(2007) study would be facilitated by discriminating between blue (star-forming) and red 
(non-star-forming) populations. Given the above discussion, it is likely that their result for 
the lowest mass bin is an underestimate. The disagreement seen in the lowest mass bin in 
Figure 6 is therefore not surprising. To summarize, while for the intermediate mass galaxies 
the SFRDs at different redshifts are in excellent agreement, the observed discrepancy for 
the lowest mass galaxies is likely caused by incompletness and selection bias in the Zheng et 
al sample. Correcting for these effects brings the Zheng et al's result into better agreement 
with that in the present study. For the highest mass bin (M > 10^^'^^ Mq), there is serious 
incompletness in our sample, which explains the relatively lower SFRDs measured for 
galaxies with higher stellar mass from the present study. 

We find the number density evolution in the present sample to depend on the stellar 
mass of galaxies. The impact of this number density evolution on the SFRD can be 
interpreted by considering the ratio of the SFRD (p*) to the number density of star forming 
galaxies (pn), P*/pn, defined as the "characteristic star formation rate" (cSFR), as a 
function of mass and redshift. The cSFR remains flat for our low-mass galaxies, at a level 
of ~ 1 Mq yr~^, while dropping by more than an order of magnitude for high mass systems 
from ^ 2M0yr^"'^ at z ~ 1 to ~ O.lMQyr^"*^ at z ~ 0.3. Intermediate masses show an 
intermediate level of evolution, with their cSFRs comparable to the high-mass objects at 
high-2;, and to low-mass objects at \ow-z. This behavior exhibits the characteristics of 
"downsizing" in galaxy evolution in the sense that the cSFR decreases with redshift faster 
for massive galaxies than for low-mass systems. 

Using spectroscopic observation s of galaxies o v er the range 1 < z < 2, selected from 



Gemini Deep Deep Survey (GDDS), 



Juneau et al. 



fl2005f ) studied changes in the SFRDs 
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with redshift as a function of stellar mass of galaxies. This provides a natural extension of 
the COSMOS study to 2; ~ 2. To study the behavior of this relation over the redshift range 
covered by the c ombined COSMOS and GDDS, we divide our sample into the same mass 
bins as those in 



Juneau et al. 



(120051 ) and compare the results in Figure [71 For the most 
massive galaxies we find the trend in SFRD with redshift to continue to z ~ 2, however 
for the less massive systems it flattens around z ~ 1.1. This implies that downsizing has 
been effective at 2; > 1, and that the massive galaxies were formed before z ~ 2 after going 
through a period of intensive star formation activity. Results in Figures 5 and 7 confirm 
that, at any given redshift out to ^ ~ 2, massive galaxies contribute less to the total SFRDs 
than objects with smaller stellar mass. However, although the sample here is selected to be 
close to Juneau et al's galaxies (in terms of the range in their stellar mass), differences in 
the selection criteria between the two samples is likely to affect the above result. 



5.3. Evolution of the SFRD as a Function of Stellar Mass and Galaxy Type 

The coverage of the COSMOS field by the Advanced Camera for Surveys (ACS), 
on-board the Hubble Space Telescope (HST), provides high resolution imaging of galaxies 
and measurement of their morphologies. Since the coverage is in only one HST-ACS band, 
it is not possible to measure morphologies for galaxies at the same rest-frame wavelength 
in all the 4 redshift bins studied here. We therefore use the spectral types as a proxy 
to rest-frame ACS morphologies. We have already shown that our estimated spectral 
types closely agree with the rest-frame morphologies in COSMOS (Capak et al 2007b). 
Nevertheless, one needs to be cautious that the spectral types and morphologies of galaxies 
are only loosely related. In the following discussion, these terms are used synonymously. 
Moreover, relations involving photometric redshifts, spectral types and stellar masses of 
galaxies should be interpreted carefully as these parameters are not independently derived 
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here. 

In Figure [HI we present the evolution of the SFRD with redshift as a function of both 
stellar mass and spectral type of galaxies. The spectral types are divided into those with 
redder SEDs (early types), intermediare color SEDs (spirals) and bluer SEDs (late- types 
and starbursts)- (Mobasher et al 2007). As with any type classification, there will be some 
level of misclassification or systems that do not fall easily into the defined spectral types. 
For the following discussion this limitation must be kept in mind. In particular, there is 
some ambiguity regarding irregular galaxies, which may in some instances fall either in 
the "spiral" or "starburst" classes as all these are represented by SEDs similar to those 
of star-forming galaxies. Nevertheless, we find a clear trend for all the spectral types in 
Figure 8. Due to serious incompleteness in galaxies with redder SEDs at log(M/M0) ~ 9.5, 
the relation for these galaxies is not shown in Figure 8. For the low and intermediate 
mass galaxies, starbursts dominate the star formation activity at high redshifts with a 
comparable contribution from spirals at lower redshifts. Also, for the lowest mass systems 
(9.5 < log{M / Mq) < 10), we find a steep increase with redshift in the SFRD for starbursts, 
compared to that for the spirals which remain unchanged with redshift. However, for the 
intermediate mass systems (10 < log{M / Mq) < 10.5), the slope of the SFRD-z relation 
increases for both spirals and starbursts, with the spirals making a relatively higher 
contribution to the total SFRD at lower redshifts. For higher mass galaxies, the sample 
is dominated by early-types, with the spirals showing a relatively steeper change in their 
SFRD with redshift compared to earlier types. This is consistent with the expectation from 
the downsizing scenario that the most massive systems are dominated by early-type old 
galaxies, which have undergone intense star formation activity at higher redshifts (z ~ 1). 
An interesting feature in Figure 7 is that the slope of the SFRD-z relation for different 
spectral types of galaxies are strong functions of their stellar mass. 
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Perhaps the most reveahng population in this diagram are the starbursts. These are a 
dominant form of star formation in galaxies at z > 1, but their prevalance at lower redshifts 
declines progressively with decreasing redshift. This is reflected in the sharp decrease from 
2; ~ 1 to 2 ~ 0.20 of the starburst contribution to the SFRD with redshift for low and 
intermediate mass systems. The fact that most definitions of a starburst favor low-mass 
star forming galaxies can be seen in the relatively large contribution of this population in 
the lowest-mass bin, and progressively less towards higher masses. In the high-mass bin, 
there are too few starbursts observed to reliably fit luminosity functions, and the lower 
limits shown refiect the sum over the detected population only. Overall, starbursts account 
for the majority of the star formation in low- to mid-mass systems at higher redshifts 
{z > 0.7), but at lower redshifts the dominant population shifts to the spiral types in these 
mass ranges. The lower the galaxy mass, the lower the redshift to which the starbursts 
remain dominant, consistent with the low- mass irregular starburst galaxies seen in the local 
universe (Figure 3 in Kennicutt 1998). This is a consequence of the change of slope of 
SFRD- 2; relations found in different stellar mass intervals. 

The spiral types are the dominant contributors at low redshift and low and intermediate 
mass, although they make a contribution similar to the early types in the highest mass bin. 
The evolution of the SFRD for spirals is flatter than any other classes, which is due to a 
transformation from starburst to spiral spectral types with decreasing redshift. We should 
note that, particularly for early types, the lowest mass systems are expected to form their 
stars very late in most hierarchical models. 



6. Discussion 
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as a function of the stellar mass of galaxies (iGlazebrook et al 
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Zheng et al.ll2007l ). Here we investigate this relation using a large and homogeneous sample 
of galaxies with well-known selection biases, accurate photometric redshifts, SFR estimates 
and stellar mass measurements. 

We find a strong dependence of the SFRD on stellar mass in the range 0.20 < z < 1.10, 
with higher mass galaxies (M > lO^*^'^ M©) contributing less (by an order of magnitude) 
to the global SFRD than lower mass systems at any given redshift. However, the rate of 
evolution of the SFRD with redshift is the same, regardless of the stellar mass of galaxies, 
in agreement with the independent study by Zheng et al (2007), in which the effect of dust 
obscuration in selection of star forming galaxies and measurement of their SFRs is taken 
into account. 



Using the GDDS with spectroscopic redshifts in the range 1 < z < 2, lJuneau et al. 



(120051 ) studied this relation from ^ ~ 1 to 2; ~ 2. The combined COSMOS and GDDS data 
show that for high mass galaxies (11.06 < log{M / Mq) < 11.76), the upward trend in the 
SFRD continues to 2 ~ 2, while for lower mass systems, after a steep increase with redshift, 
the SFRD flattens at z > 1. This implies that massive galaxies produced most of their stars 
by 2; ~ 2 (when the Universe was ~ 3.5 Gyr old), with the less massive galaxies hosting 
efficient star formation activity only after 2; ~ 1. These results confirm that the stellar mass 
of galaxies regulates the relative contribution of a galaxy to the global SFRD, consistent 
with the "downsizing" scenario for the formation of galaxies. Moreover, we find that galaxy 
"downsizing" was already occuring at 2; ~ 2 and has continued to the present epoch. 

The estimated SFRD in this study is likely to be a lower limit due to selection at optical 
bands biasing against dusty star-forming galaxies, which could significantly contribute to 
the global SFRD. Using a near-IR selected sample with Kg < 21.5, mimicking a mass 
selected sample complete to M ~ 10^° M0, Caputi et al (2006) estimated the space density 
of Luminous Infrared Galaxies (LIRGs) and Ultra-Luminous Infrared Galaxies (ULIRGs) 
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in the range 0.5 < z < 1, similar to the redshift range covered in the present study. They 
found that 24% of galaxies with lO^M© < M < 2.5 x 10" M© are either LIRGs or ULIRGs 
(excluding AGNs), with a negligible fraction of LIRGs and ULIRGs to have masses in 
excess of 2.5 x 10^^ Mq. Furthermore, using the MIPS observations in the GOODS-S field, 
Caputi et al (2005) show a very small fraction of LIRGs or ULIRGs at 0.4 < z < 1 to have 
M < 5 X 10^° M(T^. Indeed , ULIRGs of any stellar mass are found t o be very rare below 



z ~ 0.5 flFlores et al 



1999 



Caputi et al 



20061: 



Le Floc'h et al. 



20051 ) . The conclusion from 



these studies is that, the fraction of massive (M > 10^^ Mq) LIRGs and ULIRGs at 2; < 1 
is very small, implying that these objects do not significantly contribute to the total SFRD 
over this redshift range. Since the majority of the LIRGs and ULIRGs at 0.5 < z < 1 are 
associated with intermediate mass systems, the absence of these galaxies from our sample is 
likely to lead to an underestimation of the cont ribution from these ob jects (i.e. intermediate 



Le Floc'h et al 



20051) found that at z > 0.6, 



mass systems) to the global SFRD. Moreover, 
dusty starbursts with SFR > 10 Mq yr~^ dominate the cosmic star formation density, with 
their contribution increasing from 50% at z ~ to 80% at 2; ~ 1 (i.e. only 20% of the 
SFRD at 2; ~ 1 is due to the UV bright galaxies). Any selection bias against these galaxies 
may possibly compromise our estimate of the SFRD in the range 0.6 < z < 1.2. We believe 
that a possible bias against such galaxies does not affect our results here. By estimating 
the SFRDs through integrating the rest-frame UV(2800 A) luminosity functions (corrected 
for dust obscuration) over their entire magnitude range, we incorporate a contribution from 
galaxies at the faint-end of the UV luminosity function, where heavily obscured starbursts 
would lie. The excellent agreement between our estimated SFRDs from COSMOS and those 
based on the sample compiled by Hopkins and Beacom (2006) in Figure 4 supports the fact 
that the bias against dusty starburst galaxies does not significantly affect the results in the 
present study. We note that Hopkins and Beacom (2006) estimate the total contribution 
to the SFRD from their compilation by combining the SFRD inferred from the UV and 
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far-infrared luminosities. Nevertheless, the limitations imposed by selection criteria and any 
biases so introduced remains a major source of uncertainty in any such study. 

Using a sample of 24 //m selected galaxies in GOODS-S, Daddi et al (2007) found a 
roughly linear relation between the stellar mass and SFR densities. For a given mass, the 
SFRD at ^ ~ 2 was larger by factors of ~ 4 and ~ 30 relative to the star formaing galaxies 
at z ~ 1 and z ~ respectively. This is roughly consistent with our results in the two 
highest mass bins in Figure 6. In the range 1 < z < 2, the observed 24 //m wavelength 
corresponds to rest-frame 8 — 12 /xm bands, sampling the PAH features which are most 
sensitive to the star forming galaxies. The agreement between or results and Daddi et al 
(2007) further implies that the potential selection bias against dusty star forming galaxies 
in our sample is not significantly affecting the estimated SFR densities. 

By dividing our sample into three spectral types; ellipticals, spirals and starbursts, we 
study the evolution (with redshift) of the SFRD-Mass relation for each type. A striking 
feature, presented in Figure 7, is that the slope of the SFRD— ^ relations, found for different 
spectral types of galaxies are strong functions of the stellar mass. We find the massive 
galaxies, which are dominated by early-types, to significantly contribute to the total SFRD 
at earlier epochs (2; ~ 1). This is expected as more massive galaxies are able to support 
higher rates of star formation than lower mass systems, all else being equal. For the 
intermediate mass systems, the contribution from early-type galaxies to the global SFRD 
is minimal, with no significant change with redshift. The relation for lower mass systems, 
which are dominated by spirals and starbursts, show a progressive build up of the stellar 
mass (i.e. increasing SFRD) with cosmic time. The results here imply that the stellar mass 
of galaxies plays a crucial role in determining their star formation activity and spectral 
types and governs the contribution of individual galaxies to the total SFRD. 
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7. Summary and Conclusion 

Using multiwavelength data from the COSMOS survey we studied the dependence of 
the cosmic star formation density on redshift, stellar mass and spectral types. We explored 
the main parameters which govern the star formation process in galaxies. Our main results 
are summarized as follows: 

• There is a strong dependence of the SFRD on stellar mass of galaxies, with the most 
massive syatems (log(M/M0) > 10.5) contributing least (by a factor of ~5) to the 
cosmic SFRD at any given redshift. 

• Combining data from the COSMOS and GDDS, we extend the SFRD — mass 
relation to 2; ~ 2. For high mass galaxies, we find a steep increase in this relation till 
z ~ 2. This implies that the massive galaxies seen today, went through intensive star 
formation activity at ^ > 1, to generate their current stellar mass. For the less massive 
systems, the SFRD — z relation flattens at 2; ~ 1, indicating that these systems are 
currently undergoing build up of their mass, with most of their stellar mass formed at 
lower redshifts {z < 1). 

• We study dependence of the SFRD— 2;— mass relation on the spectral type of galaxies. 
We find that the slope of the SFRD— 2; relation for different spectral types is a strong 
function of their stellar mass. For low and intermediate mass systems, the main 
contribution to the cosmic SFR comes from the star-forming galaxies while, for more 
massive galaxies, older and redder galaxies are the most dominant population, with 
their contribution to the global SFRD increasing with redshift. 
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Fig. 1. — Plot of the conversion factor, /C28005 between UV luminosity (at 2800 A) and 



ongoing SFR function of redshift. L 
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The star 



formation histo ry used is taken fro r n Str olger et al. (2004), while stellar synthesis models 



are taken from 



Bruzual &: Chariot 



(l2003l ). assuming solar metallicity and a Salpeter IMF. 
The conversion factor is given in units of 10^^. 
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Fig. 2. — Distribution of the stellar mass of galaxies in different redshift intervals. The 
redshift bins are adopted following Scoville et al (2007b). The mass completeness limits for 
the COSMOS are estimated and used to select a survey complete in stellar mass. 
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Table 1: SFRDs estimated from the sum of the SFRDs for individual galaxies and from 
integrating the rest-frame Luminosity functions at four redshift intervals considered here 

z log{SFRD) Me/yr/Mpc^ log{SFRD) Me/yr/Mpc^^ 

LF simple sum 



0.315 ±0.115 
0.540 ±0.110 
0.765 ±0.115 
0.990 ±0.110 



-1.4286 ±0.0058 
-1.2733 ±0.0045 
-1.0192 ±0.0029 
-0.9656 ±0.0026 



-1.5550 ±0.0058 
-1.4196 ±0.0045 
-1.2273 ±0.0029 
-1.1033 ±0.0026 
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Fig. 3. — Distribution of the estimated SFRs (corrected for extinction) for COSMOS galaxies 
in different redshift intervals. 
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Fig. 4. — The extinction-corrected SFR densities, SFRD2800, vs. redshift, derived from COS- 
MOS galaxies are compared with other independent estimates of the SFRD. Total SFRDs 
are estimated from integrating the rest-frame UV luminosity function (filled circles) or are 
estimated by summing over the SFRs for individual galaxies (Open circles). Vertical error- 
bars correspond to Poisson counting statistics, and are all smaller than the symbol size for 
the COSMOS galaxies (blue symbols). The filled black symbols are the compilation of data 
from literature (Hopkins & Beacom 2006), reduced to the same cosmology used here. 
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Fig. 5. — Extinction corrected SFR28Q0 density from COSMOS is plotted against redshift in 
different stellar mass intervals. Errorbars are again estimated assuming Poisson statistics, 
and typically smaller than the symbol size. 
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Fig. 6. — The SFRDs from the COSMOS (this study)- (filled dots) are compared with those 
from Zheng et al (2007)- (stars). Different colors indicate different mass intervals. The 
two samples are binned in the same stellar mass intervals. The SFRDs from Zheng et al 
(2007) are estimated using the combined UV and total infrared (8-1000 /im) flux of galaxies. 
Errorbars correspond to poisson statistics. 
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Fig. 7. — Changes in SFRD w ith redshift in different stellar mass intervals, extended to 2; ~ 2 



using the 



Juneau et al 



20051) estimates from GDDS (ast erix). The SFRDs fr om COSMOS 



are re-binned to the same mass intervals as those in the I Juneau et al. 



(120051 ) sample. The 



data from GDDS is reduced to the same cosmology and IMF as for the COSMOS. 
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Fig. 8. — SFRD as a function of stellar mass and galaxy type. Blue: Starburst; Green: 
Spiral (Sa-Sd); Red: Early (E/SO); Black: all types. 



